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INTRODUCTION 

The  movement  of  water  in  soli  systems  has  long  been  of  interest  for 
applications  such  as  road,  dam,  and  reservoir  design;  agronomy;  and  ground 
thermal  storage.  Recent  proposals  for  buried,  chilled  pipelines  raise  the 
possibility  of  water  movement  in  freezing  soils  over  very  long  periods  of 
time.  The  need  to  apply  soil  water  and  ice  relations  to  practical  prob¬ 
lems  such  as  frost  heave  and  thaw  weakening  has  prompted  an  expansion  of 
research,  since  rules  of  thumb  have  repeatedly  failed.  The  physics  of  the 
processes  must  be  studied  in  increasing  depth  to  rationalize  mathematical 
models . 

The  relation  of  theoretical  microscopic  soil  water  physics  to  observed 
macroscopic  parameters  has  been  evaluated  here,  emphasizing  frost  heave 
phenomena.  Experimental  work  was  carried  out  to  determine  how  water  moves 
in  isothermal  frozen  clays  and  to  measure  simultaneous  pressures  and  temp¬ 
eratures  within  the  freezing  zone  of  silts. 


SOIL  WATER  AND  ICE 

Taber  [1]  showed  that  the  major  cause  of  volume  increase  when  soil 
freezes  is  the  formation  of  segregated  ice:  ice  formed  from  water  in 
excess  of  chat  originally  in  the  soil.  A  complete,  quantitative  explana¬ 
tion  for  the  movement  of  excess  soil  water  in  response  to  a  freezing  condi¬ 
tion  is  still  sought.  Beskow  [2]  demonstrated  that  a  pressure  gradient,  or 
a  gradient  of  free  energy,  developed  within  the  water  of  fine-grained  soil 
during  freezing.  Soil  water  will  move  under  appropriate  potential 
gradients  and  an  understanding  of  these  fluxes  in  frozen  soils  is  vital  to 
explain  frost  heaving. 


ELECTE 
JUL2  2  1982 


DISTRIBUTION  STATEMENT  A 


Approved  foe  pubha  release 
Distribution  Unlimited 


24 


LON ARP INI ,  BERG,  JENKINS,  MCGAW,  NAKANO,  OLIPHANT,  O'NEILL  AND  TICE 


Most  descriptions  view  the  soli  as  a  macroscopic  continuum  with 
reference  made  to  microscopic  phenomena  at  the  length  scale  of  a  single 
pore.  Thus  basic  laws  are  cast  In  terms  of  local  averages  which  allow 
equations  to  be  written  In  continuous  variables  despite  the  microscopic 
discontinuities.  Local  microscopic  mechanical  and  thermodynamical 
equilibria  are  usually  assumed.  Macroscopic  nonequilibria,  such  as  temp-* 
erature  variation,  may  be  viewed  as  driving  potentials,  while  equilibrium 
relations  are  applied  locally.  This  approach  Is  valid  in  soil  physics 
since  events  occur  slowly  and  phases  are  In  contact  over  large  specific 
surfaces.  Soil  hysteresis  can  greatly  complicate  computations,  but  no 
frost  heave  theories  as  yet  include  it.  Latent  heat  effects  are  large  and 
may  be  computationally  troublesome,  especially  as  they  enter  via 
non-linear 1 ties ;  references  [3,  4]  give  useful  reviews.  If  phase  transi¬ 
tion  la  concentrated  within  a  macroscopic  discontinuity,  then  a  variety  of 
methods  recently  devised  for  Stefan-type  problems  may  be  used  [5-7].  If 
the  scale  of  analysis  entails  a  spatially  extended  phase  change  zone,  a 
non-step  freezing  zone  representation  follows  naturally  [8,  9], 


Conservation  principles  of  momentum,  energy,  and  mass  enter  into  a 
complete  frost  heave  model.  Experimental  evidence  shows  that  an  increase 
in  overburden  pressure  diminishes  the  heave  produced  by  forces  acting  In 
the  freezing  zone  [10];  overburden  pressure  is  defined  macroscopically  as 
the  net  force  per  unit  area  over  all  phases  and  constituents.  The  pressure 
in  any  phase  will  normally  differ  from  the  overburden  pressure,  but  may 
equal  it  in  special  circumstances.  If  the  pore  water  gage  pressure  in  a 
saturated,  coarse-grained,  unfrozen  soil  is  zero,  then  a  weight  upon  the 
soil  must  be  supported  entirely  by  an  ’’effective  stress”  transmitted 
through  mineral  intergranular  contacts.  Microscopically  the  effective 
stresses  are  complex  and  vary  abruptly,  but  macroscopically  they  manifest 
themselves  simply  as  a  net  intergranular  force  per  unit  area  of  soil. 


The  overburden  pressure  is  not  balanced  by  an  effective  stress 
reaction  alone  but  also  by  the  pore  content  stress.  The  macroscopic  parti¬ 
tioning  of  pore  stress,  or  ”neutral  stress,"  by  the  pore  constituents 
results  from  a  microscopic  interplay  of  constituent  stresses,  interfacial 
energies,  and  geometry.  Miller's  [11]  suggestion  that  the  use  of  a  stress 
partitioning  factor  to  quantify  the  relative  participation  of  macroscopic 
pore  constituent  stresses  is  experimentally  supported  by  [12].  Using 
thermodynamics,  equations  were  derived  relating  overburden  pressure,  ice 
pressure,  liquid  water  pressure  and  volume  ratios  of  ice  and  liquid  [13, 
14].  Agreement,  however,  has  not  yet  been  reached  on  the  quantitative  dis¬ 
tribution  of  stress  within  the  phase  change  regions. 

One  may  question  the  validity  of  viewing  stresses  as  having  pressure- 
like  manifestations  macroscopically,  especially  when  the  stresses  refer  to 
an  ensemble  of  grains,  or  act  within  extremely  thin  adsorbed  water  layers. 
It  has  been  suggested  that  film  water  can  pull  nearby  water  into  it  while 
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simultaneously  pushing  away  the  ice,  if  the  water  is  stressed  aniso tropi¬ 
cally  (IS,  16].  However,  explanations  of  the  surface  layer  forces  which 
separate  a  particle  from  nearby  ice  can  be  advanced  without  recourse  to 
peculiar  water  or  stresses  [17-19].  Data  on  the  film  water  may  not  be 
necessary,  if  it  is  locally  in  equilibrium  with  the  experimentally 
accessible  pore  water.  Osmotic  pressure  also  enters  into  the  film  water  - 
pore  water  -  ice  pressure  relationships  [20].  However,  no  frost  heave 
model  takes  account  of  the  effects  of  solutes  nor  the  effects  due  to  ions 
in  clays. 

Ice  can  maintain  a  complex  state  of  stress  and  develop  stress  concen¬ 
trations,  although  it  tends  to  "relax"  and  relieve  such  conditions.  The 
idea  of  simple  ice  pressure  assumes  that  events  happen  slowly,  allowing 
micro-variations  of  stress  to  be  relieved.  Since  the  unfrozen  water  attains 
a  locally  uniform  state  and  the  water  coats  the  pore  ice  which  forms  in 
equilibrium  with  it,  the  ice  will  also  tend  to  uniformity.  When  the  soil 
matrix  and  unfrozen  water  react  to  an  imposed  load,  they  may  move  in  any 
direction,  equalizing  the  apparent  state  of  stress  in  all  constituents. 

The  pore  ice  may  also  be  considered  to  be  Isotropically  stressed,  behaving 
as  a  liquid- like  component,  with  any  macroscopic  anisotropy  entering 
through  the  stress-strain  characteristics  of  the  granular  skeleton. 


PHASE  COMPOSITION 

The  water  in  a  soil  does  not  all  freeze  at  a  unique  temperature; 
rather,  the  unfrozen  water  content  decreases  with  temperature  [21],  and 
some  temperature  drop  below  0°C  is  required  before  pore  ice  can  exist.  The 
finer  grained  a  soil  is,  the  greater  the  freezing  point  depression,  and  the 
more  unfrozen  water  it  will  retain  under  a  given  set  of  conditions. 

Williams  and  Burt  [22]  measured  the  hydraulic  conductivity  of  frozen  soil 
using  a  permeameter  with  an  aqueous  solution  of  lactose.  However,  the 
lactose  complicates  the  interpretation  of  the  data.  To  overcome  this 
difficulty  Horiguchi  and  Miller  [23]  measured  the  hydraulic  conductivity  of 
frozen  soil  in  the  temperature  range  O’C  to  -0.15®C  using  pure,  supercooled 
water.  These  experiments  quantify  the  moisture  flux  in  frozen  soils. 

The  partially  frozen  nature  of  heaving  soil  accounts  for  the  important 
unfrozen  water  flow  occurring  within  the  warmest,  least  frozen  portion  of 
the  frozen  soil.  Evidence  shows  ice  segregation  occurring  within  the 
frozen  zone,  behind  a  "frozen  fringe”  [24,  23].  This  fringe  may  be  defined 
as  the  zone  over  which  volumetric  ice  content  Increases  from  zero  at  the 
freezing  front  (0*C  isotherm)  to  100Z  at  the  warmest  lens.  High  overburden 
pressures  and  low  temperature  gradients  tend  to  increase  the  fringe 
thickness. 
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An  equilibrium  relation  which  appears  applicable  to  the  ice- liquid 
mixture  in  frozen  soil  is  the  generalized  Clapeyron  equation  [26]: 

L(T-T  ) 

<pw-»>vw  -  — r2-  •  <l) 

o 

Applied  microscopically,  this  equation  can  be  used  to  explain  equilibrium 
between  pore  ice  and  unfrozen  water.  Maeroscoplcally,  py  and  t  this 
equation  presumably  pertain  to  measured  pore  water  values.  The  Clapeyron 
equation  has  been  linked  with  another  relation  which  is  usually  presented 
in  terms  of  soil  capillarity.  When  different  water  phases  co-exist  with 
curved  phase  interfaces  (as  in  soil  pores),  there  is  a  pressure  jump 
between  the  two  phases  due  to  interfacial  effects. 


"j  -  "v  ■  v(8>  (2) 

where  j  refers  to  ice  in  saturated,  frozen  soil,  or  air  in  unsaturated, 
unfrozen  soil.  F( 9)  is  a  macroscopic  function  of  unfrozen  water  content, 
which  corresponds  to  the  mean  curvature  of  the  microscopic  phase  interface. 
Koopmans  and  Miller  [27]  hypothesized  that  F(9),  at  a  fixed  water  content, 
was  the  same  for  an  unsaturated,  unfrozen  soil  and  a  saturated,  frozen  soil. 
Pw  is  zero,  then  from  Eq.  (2) 


p. /p  -  a.  /a 
ri  ra  iw  aw 


(3) 


For  the  same  granular  soil,  Koopmans  and  Miller  [27]  found  the  pressure 
ratio  to  be  about  0.45,  indicating  that  capillary  effects  dominate  in  the 
same  way  for  both  freezing  and  drying.  For  a  clay  the  ratio  was  unity. 

When  a  particular  unfrozen,  unsaturated  soil  is  considered,  the  curves 
of  Eq.  (2)  are  often  called  soil  moisture  characteristic  curves.  Typically, 
the  curves  differ  somewhat  for  drying  or  watting  processes  (hysteresis). 

At  atmospheric  air  pressure,  neglecting  hysteresis,  one  obtains  a  curve  re¬ 
lating  moisture  content  and  liquid  pressure.  With  a  frozen,  saturated 
soil,  one  may  construct  freezing  characteristic  curves  by  relating  3  to 
Py  and  T  using  Eqs.  (1)  and  (2).  For  a  fixed  value  of  py,  a  relation 
is  found  between  3  and  T.  Thus  unfrozen  water  relations,  discussed  earli¬ 
er,-  are  special  cases  of  freezing  characteristic  curves  at  constant  pres¬ 
sure.  One  cannot  obtain  liquid  pressure  in  a  freezing  soil  directly  from 
an  unfrozen  moisture  characteristic  curve.  F(9)  may  be  the  same  in  both 
cases,  but  the  ice-water  surface  tension  Ojy  must  be  used  in  place  of 
°ay  and  the  product  aiv  F(9)  only  provides  values  for  the  difference 
between  ice  and  liquid  pressures.  The  ice  pressure  need  not  be  zero  when 
tbs  overburden  pressure  is  small  because,  in  general,  any  number  of  combin- 
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ation*  of  liquid  and  lea  prassura  nay  hold  for  given  values  of  6  and  over¬ 
burden  pressure.  Data  do  not  support  application  of  the  same  characteris¬ 
tic  curve  data  to  both  frozen  and  unfrozen  soils,  if  both  soils  are  unsatu¬ 
rated  [20].  In  the  freezing  of  unsaturated  soil,  three  pore  phases  are 
present,  and  a  more  complex  relation  between  pore  contents  must  hold  than 
that  in  Eq.  (2);  navertheless,  characteristic  curves  from  unfrozen,  unsatu¬ 
rated  soil  are  generally  used  in  frost  heave  models  of  unsaturated,  freez¬ 
ing  soil. 

The  movement  of  moisture  in  soil  systems  is  usually  described  by  some 
fora  of  Darcy's  law.  This  expresses  the  conservation  of  momentum  when 
inertial  effects  are  negligible.  The  gradient  of  free  energy,  or  thermo¬ 
dynamic  potential,  represents  a  net  force  balanced  by  the  drag  of  other 
components  on  the  liquid,  which  is  proportional  to  the  flow  rate.  Thus  the 
moisture  flux  can  be  written 

q  -  -k  7f  (4) 

The  difference  in  free  energy  of  the  soil  moisture  is  the  work  required  to 
change  the  water  from  one  state  to  another.  This  gives  the  clearest 
physical  picture  of  the  driving  potential  for  soil  water  movement  since  it 
contains  all  of  the  component  potentials  usually  encountered,  e.g. 
pressure,  gravity  electrostatic  forces,  osmosis.  The  free  energy  can  be 
expressed  as  a  function  of  external  force  fields  and  contact  forces.  Thus 
it  is  common  to  write  the  motive  potential  directly  in  terms  of  pressure 
and  gravity  fields: 


<  -  7  2)  (5) 


The  soil  suction  p  is  a  combination  of  hydrostatic  pressure,  surface 
tension  and  adsorption  forces  of  the  soil  skeleton.  The  soil  suction  p  can 
be  expressed  as  a  function  of  water  content  and  temperature,  as  noted 
already.  Then  the  extended  Darcy's  law  is 


f-  -  -De  7  8  -  D-  7T  -  k  e  (6) 

pw  1 

Eq.  (6)  may  be  substituted  into  the  conservation  of  mass  equation  to  arrive 
at 

|£  -  7[De70  +  dt7T  +  k  aj  (7) 

No  frost  heave  models  have  used  the  mass  flux  relation  as  a  function  of 
both  pressure  and  temperature  gradients,  exemplified  by  Eq.  (6).  If  one- 
dimensional  flow  is  considered,  with  negligible  temperature  effects,  then 
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In  a  number  of  frost  heave  models,  soil  moisture  dlffuslvity  is  used 
in  the  seme  manner  for  frozen  and  unfrozen  soli,  evidently  by  analogy 
between  drying  of  unsaturated  soil  and  freezing  of  saturated  soil.  How¬ 
ever,  the  analogy  may  not  hold.  In  fact,  the  Df  term  can  be  dominant  in 
saturated,  granular  soil.  Unlike  pa,  the  ice  pressure  will  not  be 
constant  in  space,  but  should  vary  from  some  small  value  where  ice  content 
is  low  to  the  overburden  pressure  where  segregation  is  occurring.  According 
to  Eq.  (2),  the  liquid  pressure  gradient  is  then 

3pw_3pi  „  rr3Fi  36  *  r3Fi  31 1 

IT  TT  aiwt  Wt  +  We  TJJ  (9) 

Thus  the  liquid  pressure  gradient  depends  upon  ice  pressure  gradient,  and 
gradients  in  0  and  T.  No  simple  relation  for  dlffuslvity  exists  unless  the 
ice  pressure  and  thermal  gradients  are  negligible. 

Pressure  or  temperature  gradients  may  cause  relative  motion  between 
ice  and  small  particles  completely  enveloped  in  the  ice  [18,  28].  The 
processes  involved  are  not  completely  understood  but  the  ice  must  melt  on 
one  side  of  a  particle,  with  liquid  being  transported  around  the  particle 
to  a  location  where  it  refreezes.  The  process,  called  regelation,  may 
operate  in  heaving  soils,  so  that  movement  of  ice  as  well  as  unfrozen  water 
may  occur  [11].  The  extent  of  such  activity  in  soil  has  not  been  estab¬ 
lished,  but  the  possibility  of  pore  ice  movement  exists. 


ISOTHERMAL  TRANSPORT  OF  WATER  IN  SUBFREEZING  SOIL  WITHOUT  ICE 

Frozen  soils  need  not  contain  ice  if  the  moisture  content  is  less  than 
the  equilibrium  value  for  the  temperature  of  Interest.  To  examine  the 
movement  of  water  under  such  conditions,  tests  were  run  whereby  wet  and  dry 
soil  columns  were  joined  and  the  water  was  observed  to  move  from  the  wet, 
ice-free  part  of  the  frozen  soil  column  to  the  dry  region  [29]. 

Assuming  the  extended  Darcy  law  is  valid,  the  movement  qf  the  soil 
water  is  governed  by  Eq.  (8),  with  the  following  conditions: 

0  x  <  0 

o 

0 

c 


9(x,Q) 


x  >  0 


(8a) 
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lia  0(x,t)  -  8 
x+—  0 

lim  8(x,c)  -  8 

Xi-H*  *" 


(8b) 


Nakano  [30]  shoved  that  a  waait  solution  of  Eq.  (8-8b)  converges  asymptoti¬ 
cally  to  a  similarity  solution  of  Eq.  (8),  assuming  that  the  soil  dif- 
fuslvlty  is  given  by 

D(8)  -  D  (8-8  )3  (10) 

o  c 


An  experimental  asymptote  can  be  obtained  at  some  finite  time,  chosen  so 
that  the  moisture  profile  no  longer  varies  with  respect  to  a  similarity 
parameter.  With  a  best  fit  of  the  experimental  and  theoretical  profiles 
the  values  of  DQ  and  8  can  be  found.  This  method  is  also  applicable  to 
the  determination  of  the  diffusivity  for  unfrozen  porous  media  and  is  a 
useful  new  technique. 

The  soil  was  Morin  clay  with  a  specific  surface  area  of  60  m2/g-soil. 
The  measured  unfrozen  water  content  (with  21Z  total  water  content)  versus 
temperature  showed  that  about  10Z  of  the  water  remained  unfrozen  at  -L.0°C. 

Graduated  plastic  pipettes  were  used  to  enclose  soil  columns  0.80  cm 
in  diameter  by  20.33  cm  in  length.  Two  soil  columns  with  the  same  soil  dry 
density  of  1.3  g/cm3  were  required  per  test,  one  uniformly  dry  and  the 
other  uniformly  wet  at  8.0Z,  The  thoroughly  mixed  water  and  soil  was 
allowed  to  set  for  at  least  three  days  to  attain  moisture  uniformity. 

An  experiment  began  by  joining  pipettes  with  dry  and  wet  soil,  both 
held  at  -1.0*C.  The  surfaces  of  contact  were  smoothed  to  avoid  contact 
resistance.  The  joined  pipettes  were  then  placed  in  a  constant  temperature 
bath  containing  an  ethylene  glycol/water  mixture  held  at  -1.0*C  and 
maintained  to  within  ±0.Q3aC.  After  a  specified  time,  the  compound  pipette 
was  removed  from  the  constant  temperature  bath  and  quickly  sectioned  into 
64  equal  segments  for  water  content  determination. 

Tests  were  carried  out  with  time  durations  of  2,  4,  8  and  16  days  and 
also  at  4,  8,  16  and  32  days,  to  examine  the  reproducibility  of  the  data. 
The  water  content  is  plotted  versus  the  distance  in  Fig.  1.  The  two  sets 
of  profiles,  after  8  days,  indicate  the  difficulty  in  reproducing  the  same 
profile.  There  was  about  0.SZ  water  content  difference  in  the  wetter  end 
of  a  soil  column.  The  water  content  8  is  plotted  versus  a  similarity  para¬ 
meter  €  (Fig.  2),  which  Indicates  that  the  experimental  profiles  tend  to 
stabilize  as  if  they  were  converging  upon  an  asymptotic  profile  reasonably 
well  after  8  days.  Using  the  least-squares  method  the  parameters  in  Eq. 
(10)  were  determined  as 
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Figure  1.  Water  movement  in  ice-  Figure  2.  Similarity  correlation  of  water 
free  Morin  clay.  9Q  »  8.0Z;  t  -  content  in  ice-free  Morin  clay.  90  «  8.0Z; 
-1.0*C.  T  -  -i.0#C. 
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D<0)  -  1.21  <9  -  .0038)*  (11) 

where  0(9)  is  in  cm2/day  and  9  is  in  g  water/g  dry  soil.  The  curve  in 
Fig.  2  is  the  computed  9  profile  using  Eq.  (11)  and  is  a  reasonably  good 
approximation,  but  is  not  quite  a  close  fit.  If  0(9)  is  evaluated  directly 
from  the  experimental  profile  at  8  days,  using  the  same  mathematical  basis, 
it  appears  that  there  is  a  slight  peak  in  D(8)  near  9  *  IT.  The  validity 
of  this  calculation  is  uncertain,  due  to  the  error  in  the  numerical  evalua¬ 
tion  of  derivatives,  but  it  is  possible  that  the  fit  is  somewhat  erroneous 
because  the  actual  0(8)  may  not  be  a  monotonically  increasing  function  of 
8. 

ISOTHERMAL  TRANSPORT  OF  WATER  IN  SOIL  WITH  ICE 

Tests  were  run  using  the  same  technique  and  apparatus  as  described  in 
the  previous  section  except  that  the  water  concent  was  adjusted  so  that  ice 
was  present  in  the  wet  section  [31].  The  total  water  content  was  20. OX  and 
the  dry  density  of  the  clay  was  1.42  Mg/m3.  The  unfrozen  water  content  was 
evaluated  by  a  pulsed  nuclear  magnetic  resonance  (NMR)  technique  to  yield 
the  following  unfrozen  water  relation  for  a  warming  soil: 
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-.5080 

9  -  .1269  (-T)  (12) 


Since  Che  tests  were  run  at  -laC,  the  initial  water  and  ice  contents  were 
12. 7Z  and  7.3Z  respectively.  Some  test  results  are  shown  in  Fig.  3,  which 
indicate  that  water  moved  slowly  but  steadily  from  the  wet  part  of  the  soil 
column- to  the  dry  part. 

Supplemental  tests  were  carried  out  to  determine  if  the  phase  equi¬ 
librium  of  the  unfrozen  water  was  maintained.  The  total  water  content  of 
the  soil  was  measured  gravimetrically  while  Che  unfrozen  water  content  was 
determined  by  NMR.  As  the  total  water  content  decreased  the  unfrozen  water 
concent  remained  at  the  initial  value  until  the  total  water  content  dropped 
to  about  14Z.  The  unfrozen  water  content  of  the  wet  column  remained 
constant  at  the  equilibrium  unfrozen  water  content  at  -1°C  and  any  net 
water  removal  from  the  wet  soil  column  was  at  the  expense  of  the  ice.  The 
ice  which  melted  in  Che  vet  column  absorbed  energy  from  the  surroundings 
but  the  system  temperature  did  not  vary  more  chan  t.01*C.  Thus  the  ice 
melted  so  slowly  chat  latent  energy  could  be  transferred  with  a  very  minor 
temperature  gradient. 

If  Darcy's  law  is  valid  for  water  flow  with  ice  Chen 
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Using  Eq.  (13)  it  was  possible  to  calculate  the  diffusivity  using  the  data 
for  9t  and  q  as 
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D(8t)  -  49.9  (8to  -  8t)*678  (14) 

By  solving  Che  water  concent  equation  with  a  dlffusivicy  given  by  Eq. 
(14),  the  cumulative  water  flux  was  predicted  as  a  function  of  time  and 
position.  Since  these  values  compared  quite  well  to  the  data,  a  Darcy  type 
law  may  govern  the  flow  of  film  water  in  icy  soils  under  conditions  similar 
to  those  of  the  experiment. 

FREEZING  ZONE  TEMPERATURE  AND  PRESSURE 

Although  moisture  tensions  are  known  to  exist  below  an  advancing 
freezing  front,  and  are  suspected  to  exist  within  the  freezing  zone  itself, 
the  cause  of  these  tensions  has  not  been  adequately  explored.  Critical 
assumptions  in  frost  heave  theories  are  not  adequately  supported  by  data  on 
temperatures  and  pressures  that  exist  within  the  freezing  zone.  Tests  were 
designed  to  provide  concurrent  measurements  of  temperature  and  pressure  in 
the  moisture  phase  of  freezing  silt  soils.  Graves  sandy  silt  and  Northwest 
silt  were  frozen  unidirectionally  from  the  upper  surface  with  the  lateral 
surfaces  insulated.  Selected  physical  properties  for  these  soils  are 
listed  in  reference  [31]. 

A  special  apparatus  was  designed  and  fabricated  to  contain  the  300-mm- 
long  soil  specimen  during  freezing  and  to  support  the  measuring  devices  for 
temperature  and  moisture  tension  [31 J.  Top  and  bottom  temperatures  were 
maintained  to  a  tolerance  of  0.01°C  by  separate  fluid  baths;  the  upper  tem¬ 
perature  was  adjusted  daily  to  advance  the  freezing  front  into  the  specimen 
at  a  generally  decreasing  rate  between  20  and  5  mm/ day  over  a  24-day  peri¬ 
od.  De-alred  water  was  supplied  continuously  to  the  bottom  of  the  soil 
through  a  porous  stone  connected  to  an  external  4*C  water  supply  held  level 
with  the  bottom  of  the  specimen. 

Temperatures  along  the  length  of  the  soil  column  were  measured  using 
calibrated  thermistors  mounted  within  the  wall  of  the  containing  vessel  in 
direct  contact  with  the*  soil.  These  registered  soil  temperatures  within 
0.05*C.  Moisture  tensions  were  measured  in  the  unfrozen  and  frozen  por¬ 
tions  of  the  soil  column  by  means  of  special  tensiometers  designed  at  CRREL 
for  temperatures  above  and  below  0*C  [32].  A  IGZ  solution  of  ethylene  gly¬ 
col  and  water  was  used  within  the  porous  cup  which  remains  unfrozen  to  a 
temperature  of  approximately  -4*C  and  does  not  diffuse  into  the  soil 
water.  Within  experimental  error  the  soil  moisture  tension  is  found  to  be 
the  same  for  the  10Z  solution  as  for  water. 

Figure  4  is  an  example  of  the  data  obtained  for  the  Graves  sandy 
silt.  The  Northwest  soil,  a  true  silt,  developed  similar  features.  Ten¬ 
sion  was  registered  continuously  by  the  tensiometers  and  differed  signifi¬ 
cantly  from  atmospheric  pressure  only  when  the  soil  temperature  was  below 


Depth  (mm) 
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(*P)  (Atm) 


-0.4*C.  The  value  of  moisture  tension  increased  rapidly  as  cemperatures 
became  lower,  finally  reaching  an  apparent  maximum  at  approximately  70  kPa, 
after  which  it  remained  constant  or  fell  off  slightly  at  a  temperature  of 
-l.5*C  for  the  two  soils  tested.  Water  flowed  through  the  unfrozen  soil  at 
a  barely  detectable  moisture  tension  gradient,  consistent  with  the  measured 
water  intake  and  the  hydraulic  conductivity  of  the  soils  tested.  Moisture 
tension  gradients  in  the  freezing  zone  reached  a  maximum  of  approximately  2 
kPa/mm. 

Figure  5a  is  a  graph  of  the  moisture  tensions  and  temperatures  mea- 
sured  in  the  freezing  zone  of  the  Graves  sandy  silt  plotted  for  various 
depths.  Figure  5b  is  a  graph  of  unfrozen  water  content  for  this  soil,  as  a 
function  of  below-freezing  temperature,  measured  with  the  NMR  procedure. 

Combining  the  data  of  Fig.  5a  and  5b  leads  to  an  important  new  step  in 
the  assessment  of  frost  heaving:  a  moisture-tension  characteristic  for  the 
unfrozen  water  behind  the  freezing  front  (Fig.  6).  The  combining  of  pres¬ 
sure  and  temperature  data  assumes  that  the  pressure  measured  by  the  tensi¬ 
ometers  is  that  of  the  unfrozen  water  itself,  i.e.  the  state  of  the  water 
is  the  same  in  the  active  freezing  zone  as  it  is  in  the  unfrozen  water 
tests.  An  observed  correspondence  in  the  two  types  of  test  provides  strong 
support  for  this  assumption:  the  temperature  of  nucleatlon  found  by  ice 
seeding  in  the  unfrozen  water  tests  (approximately  -0.5®C)  is  the  same  tem¬ 
perature  at  which  moisture  tension  first  exceeds  the  hydrostatic  level  in 


// 
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Figure  5.  Moisture  tension  (a)  and  unfrozen  Figure  6.  Freezing-zone  moisture  char- 

water  content  (b)  versus  temperature  in  freezing  acterlstic;  Graves  sandy  slit, 

zone  of  Graves  sandy  silt. 
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the  advancing  freezing  front.  Moreover,  a  steep  gradient  of  tension  devel¬ 
ops  at  temperatures  immediately  lower  than  -0.5aC. 


Unlike  a  moisture  characteristic  for  an  unfrozen  soil,  the  freezing- 
zone  moisture  characteristic  is  not  independent  of  temperature;  each  inter¬ 
section  of  moisture  tension  and  unfrozen  water  content  on  such  a  plot 
corresponds  to  a  narow  range  of  temperatures  below  0SC.  The  range  is  most 
likely  the  result  of  local  variations  in  gradation,  density,  or  water 
content  in  the  vicinity  of  a  tensiometer.  A  normal  moisture  characteristic 
corresponds  to  an  air/ water  interfacial  system,  while  the  freezing-zone 
moisture  characteristic  apparently  conforms  to  a  water/ ice  interfacial 
system  having  a  lower  surface  energy. 


CONCLUSIONS 

New  relations  for  the  flow  of  water  in  partially  frozen  soils  are  now 
being  incorporated  into  frost  heave  models.  The  thermodynamic  equilibrium 
equations  for  freezing  soil  water  are  well  established  but  the  partition  of 
the  soil  system  stresses  within  freezing  zones  is  poorly  understood. 

It  has  been  demonstrated  chat  the  hydraulic  conductivity  of  water  in 
frozen  soils  is  small  but  non-zero.  Over  time  spans  of  several  years, 
which  are  appropriate  for  new  large  scale  engineering  projects  in  cold 
climates,  water  flow  in  frozen  soils  will  be  significant.  The  experimental 
work  has  shown  that  Darcy's  law  is  valid  for  film  water  flow  in  ice-free 
frozen  soils  and  is  a  reasonable  approximation  for  frozen  soil  with  ice  if 
the  ice  pressure  gradient  is  negligible. 

The  capability  of  deriving  the  freezing-zone  moisture  characteristic 
for  a  frost-susceptible  soil,  demonstrated  here,  is  important  in  the 
theoretical  prediction  of  the  amount  of  segregated  ice  in  a  freezing  soil.  w. 
The  empirical  correspondence  between  temperature,  moisture  tension,  and  ^ 
unfrozen  water  content  may  lead  to  physically  realistic  computer  algorithms 
for  the  freezing  zone.  Quantitative  predictions  from  frost  heave  models 
will  rely  on  new  Information  of  this  kind. 
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NOMENCLATURE 
D  diffusivity 

Dg  k/pwg(-||-).j,  moisture  diffusivity 

D-j  temperature  diffusivity 

ez  unit  vector  in  z-direction 
f  free  energy 
F  .  surface  tension  function 
g  gravitational  constant 

k  hydraulic  conductivity 

L  latent  heat  of  fusion 
p  pressure 
q  liquid  flux 
t  time 
T  temperature 

T0  freezing  temperature  of  bulk  water 
v  specific  volume 


x  coordinate  in  flow  direction 

z  coordinate  in  gravity  direction 

9  liquid  water  content 
90  initial  water  content 
9t  total  water  content 
5C  x  t^ 

it  osmotic  pressure 
0  density 
a  surface  tension 
metric  potential 
V  nabla  operator 
Subscripts 

a,i,w  air,  ice,  water,  respec¬ 
tively 

aw,iw  air/water,  ice/water  re¬ 
spectively 


